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Abstract. The aim of this paper is to study the relationship between 
Hamiltonian dynamics and constrained variational calculus. We de- 
scribe both using the notion of Lagrangian submanifolds of convenient 
symplectic manifolds and using the so-called Tulczyjew's triples. The 
results are also extended to the case of discrete dynamics and nonholo- 
nomic mechanics. Interesting applications to geometrical integration of 
Hamiltonian systems are obtained. 

1. Introduction 

One of the main notions in symplectic geometry is the concept of La- 
grangian submanifolds. This concept arises in several and different interpre- 
tations of physical, engineering and geometric phenomena. In this paper, we 
will focus our attention in their applications to Lagrangian and Hamiltonian 
dynamics of constrained systems. 

For instance, the theory of Lagrangian submanifolds gives a geometric 
and intrinsic description of Lagrangian and Hamiltonian dynamics [344 [35] . 
Moreover, it allows us to relate both formalisms using as a main tool the 
so-called Tulczyjew's triple 

T*TM TT*M — T*T*M. 

Recall that, in the above expression, au is the Tulczyjew's canonical sym- 
plectomorphism from TT*M (with its canonical symplectic structure 
(Itojm) to T*TM (equipped now with its canonical symplectic structure 
w™)i and /3m is the canonical symplectomorphism defined by the sym- 
plectic structure ujm on T*M. The Lagrangian dynamics is "generated" 
by the Lagrangian submanifold dL(TM) of T*TM where L : TM -> M is 
the Lagrangian function, while the Hamiltonian formalism is generated by 
the Lagrangian submanifold dH(T*M) of T*T*M where H : T*M ->■ M 
is the corresponding Hamiltonian energy. The dynamics and the relation- 
ship between both formalisms are based on the central part of the Tulczy- 
jew's triple, TT*M, where the Lagrangian submanifolds «M(dL(TM)) and 
P]£(dH(T*M)) live. Of course, any submanifold of a tangent bundle auto- 
matically determines a system of implicit differential equations; in this case, 
we can apply the integrability constraint algorithm, described in § 12.21 (see 
[9j 132] for more details), to find the integrable part of the dynamics defined 
by this Lagrangian submanifold. 

This model is also valid for constrained variational calculus, determined 
by a function L : C — > R where C is a submanifold of TM, with inclusion 
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ic '■ C — > TM. In this case, we can also construct a new Lagrangian 
submanifold (see § HJ S L of T*TM as: 

E L = {|i£ T*TM \i* c V = dL}. 

Thus, we can obtain via cxm a new Lagrangian submanifold of the tangent 
bundle TT*M which completely determines the equations of motion of the 
constrained dynamics (see [E]), which are, in a regular case, of Hamiltonian 
type. In this paper, we will also prove that given an arbitrary Hamiltonian 
system we can construct a (possibly) constrained Lagrangian system that 
generates the original one. 

It is necessary to stress that the equations derived are purely variational 
and, consequently, different from the nonholonomic equations obtained by 
applying the Lagrange-D'Alembert's or Chetaev's principles (see 0E1E3] for 
more details). It is well-known that the nonholonomic equations give the 
right physical dynamics of a constrained mechanical systems [26], mainly 
related to the rolling motion. Geometrically, nonholonomic constraints are 
globally described by a submanifold C of the velocity phase space TM. 
If C is a vector subbundle of TM, we are dealing with linear constraints 
and, when C is an affine subbundle, we are in the case of affine constraints. 
Lagrange-D'Alembert's or Chetaev's principles allow us to determine the 
set of possible values of the constraint forces only from the set of admissible 
kinematic states, that is, from the constraint submanifold C determined 
by the vanishing of the nonholonomic constraints. An interesting study of 
nonholonomic systems as implicit differential equations is presented in |19j . 

On the other hand, constrained variational calculus is mainly related to 
mathematical and engineering applications, specially in control theory. In 
this paper, we will show the close relationship between classical Hamiltonian 
dynamics and constrained variational calculus. In fact, both are equivalent 
under some regularity conditions. Nevertheless, nonholonomic mechanics is 
described by a general submanifold of TT*M which is not Lagrangian; this 
fact implies the non-preservation properties of the nonholonomic flow. 

Moreover, we will study the discrete formalism, which will be also inter- 
preted in the same way as Lagrangian submanifolds of the cartesian product 
of two copies of T*M, equipped with a suitable symplectic structure. We are 
interested in finding a geometrical answer to two, in principle, alternative 
ways to derive geometric integrators for constrained Lagrangian problems: 
a pure discrete variational procedure or a symplectic numerical method for 
the associated Hamiltonian system. We will show in which cases both pro- 
cedures match. Additionally, we will derive a general algorithm to compare 
the dynamics in variational constrained and nonholonomic cases, both in 
continuous and discrete cases. 



2. Geometric preliminaries 

2.1. Lagrangian submanifolds. In this section we will introduce some 
particular constructions of Lagrangian submanifolds that are interesting for 
our purposes (see [27} 155]). 
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First, let us recall that given a finite-dimensional symplectic manifold 
(P, uj) and a submanifold N, with canonical inclusion : N ^ P, then A 
is a Lagrangian submanifold if i* N uj = and dim A = ^dim P. 

i) An interesting class of Lagrangian submanifolds, which will be useful 
111 § El is the following. Let (P, uj) be a symplectic manifold and g : 
P — > P a diffeomorphism. Denote by Graph (g) the graph of g, that 
is Graph (g) = {(x, g (x)) , z G P} C P x P, and by pri : P x P — >■ P, 
z = {0,1}, the canonical projections. Then (P x P, 0), where = 
prjf w — prQ uj, is a symplectic manifold. Let i g : Graph (g) P x P 
be the inclusion map, then 

= (pr )* (g*u-u) . 

Thus, 5 is a symplectomorphism (that is, g*uj = uj) if and only if 
Graph g is a Lagrangian submanifold of P x P. 

A distinguished symplectic manifold is the cotangent bundle T*M of any 
manifold M. If we choose local coordinates (q l ), 1 < i < n, then T*M 
has induced coordinates (q l ,Pi)- Denote by ttm '■ T*M — >• M the canonical 
projection defined by i^M^q) = Q, where e q £ T*M. Define the Liouville 
one- form or canonical one- form #a/ G K l T*M by 

((0 M )e , X) = (e , Ttt m {X)), where X G T € T*M , e G T*M. 

In local coordinates we obtain 6 m = Pi dg*. The canonical two- form ujm on 
T*M is the symplectic form ujm = —&@M (that is um = d(f A d^). 

ii) Now, we will introduce some special Lagrangian submanifolds of 
the symplectic manifold (T* M,ujm)- For instance, the image T,\ = 
X(M) C T*M of a closed one-form A G A 1 M is a Lagrangian sub- 
manifold of (T* M,ujm), since A*ojm = — dA. We then obtain a sub- 
manifold diffeomorphic to M and transverse to the fibers of T*M. 
When A is exact, that is, A = d/, where / : M — > R, we say that 
/ is a generating function of the Lagrangian submanifold Y,\ = Ej. 
Locally, this is always the case. 

A useful extension of the previous construction is the following result due 
to W.M. Tulczyjew. 

Theorem 2.1 ([S],[35]). Let M be a smooth manifold, t m : TM M its 
tangent bundle projection, N C M a submanifold, and f : N — ► M. T/ien 

E/ = {p G T*M I tt m (p) G TV and (p,u) = (d/,u) 

/or a// v G TA C TM such that tm(v) = t^m(p)} 
is a Lagrangian submanifold ofT*M. 

Taking / as the zero function we obtain the following Lagrangian sub- 
manifold 

E = {p G T*M\ N | (p, v) = 0, Vu G TA with tmW = 7r M (p)} , 
which is just the conormal bundle of A: 

^) = {eer*M| JV ;e| w = 0}. 
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Given a symplectic manifold (P,oj), dimP = 2n it is well-known that its 
tangent bundle TP is equipped with a symplectic structure denoted by dpto 
(see |25j). If we take Darboux coordinates (q l ,Pi) on P, 1 < i < n, then 
uj = dq l A dpi and, consequently, we have induced coordinates (q l ,pf, q l ,Pi), 
(q l ,Pi, di, b l ) on TP and T*P, respectively. Thus, dpu: = dq l A dpi + dq l A dpi 
and up = dq l A doj + dpi A db\ If we denote by b w : TP — )• T*P the iso- 
morphism defined by w, that is b w (t>) = i v uj, then we have \>u(q l ,Pi; q l ,pi) = 
(q l ,Pi, —pi, q 1 )- Given a function H : P — > R, and its associated Hamiltonian 
vector field that is, ix H uJp = dH, the image Xh{P) is a Lagrangian 
submanifold of (TP, dpujp). Moreover, given a vector field X € 3£(P), it 
is locally Hamiltonian if and only if its image X(P) is a Lagrangian sub- 
manifold of (TP, dpuj). It is interesting to note that dpuj = — b^wju and 
K(X H (M)) = dH(M). 

As it is briefly mentioned above, an important notion in the theory of 
Lagrangian submanifolds is the concept of generating function. If we have 
a Lagrangian submanifold N of an exact symplectic manifold (P, u) = dO), 
where 9 G A 1 P, then = i* N co = i* N d9 = d(i* N 9). Consequently, applying 
the Poincare's lemma, there exists a function S : U — > R defined on a 
open neihborhood U of N such that i* N 9 = dS. We say that S is a (local) 
generating function of the Lagrangian submanifold N . 

2.2. Implicit differential equations. An implicit differential equation on 
a general smooth manifold M is a submanifold P C TM. A solution of P 
is any curve 7 : / — > M, / C R, such that the tangent curve (7(i),7(i)) € P 
for all t € I. The implicit differential equation will be said to be integrable 
at a point if there exists a solution 7 of E such that the tangent curve passes 
through it. Furthermore, the implicit differential equation will be said to be 
integrable if it is integrable at all points. Unfortunately, integrability does 
not mean uniqueness. The integrable part of E is the subset of all integrable 
points of E. The integrability problem consists in identifying such a subset. 

Denoting the canonical projection tm '■ TM — > M, a sufficient condition 
for the integrability of E is 

E C TM, 

where C = tm(E), provided that the projection tm restricted is a submer- 
sion onto C . 

2.2.1. Extracting the integrable part of E. A recursive algorithm was pre- 
sented in [32j that allows to extract the integrable part of an implicit differ- 
ential equation E. We shall define the subsets 

Eq = E, Co = C, 

and recursively for every k > 1, 

Ek = Ek-i n TCk~i, Ck = Tu(Ek), 

then, eventually the recursive construction will stabilize in the sense that 
E k = E k+1 = ... = Eoo, and C k = C k+1 = ... = C^. It is clear by 
construction that Eoo C TCoq. Then, provided that the adequate regularity 
conditions are satisfied during the application of the algorithm, the implicit 
differential equations Poo will be integrable and it will solve the integrability 
problem 
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3. TULCZYJEW'S TRIPLES 

In this section we summarize a classical result due to W.M. Tulczyjew 
showing a natural identification of T*TM and TT*M, where M is any 
smooth manifold, as symplectic manifolds. This construction plays a key 
role in Lagrangian and Hamiltonian mechanics. 

Is easy to see that TT*M, T*TM and T*T*M are naturally double vec- 
tor bundles (see [12], [33]) over T*M and TM. In [34] and [35], Tulczyjew 
established two identifications, the first one between TT*M and T*TM (use- 
ful to describe Lagrangian mechanics) and the second one between TT*M 
and T*T*M (useful to describe Hamiltonian mechanics). The Tulczyjew 
map ajf is an isomorphism between TT*M and T*TM. Beside, it is also 
a symplectomorphism between these double vector bundles as symplectic 
manifolds, i.e. (TT*M , (Itujm), where drWM is the tangent lift of um, and 
(T*TM,ujtm)- In the following diagram we show the different relationships 
among these bundles. 



TT*M T*TM 




M 



The definition of T*tm is given in the following remark. 

Remark 3.1. Given a tangent bundle tn : TN — > N, for each y G T X N we 
can define 

V y = ker {T y T N : T y TN -> T X N} , r N (y) = x. 

Summing over all y we obtain a vector bundle V of rank n over TN. Any 
element u £ T X N determines a vertical vector at any point y in the fibre 
over x, called its vertical lift to y, denoted by u v (y). It is the tangent vector 
at t = to the curve y + tu. If X is a vector field on N, we may define 
its vertical lift as X v (y) = (X(r N (y))) V . Locally, if X = X^ in a local 
neighborhood U with local coordinates x % , then X v is locally given by 

X y = x^, 

ov l 

with respect to induced coordinates {x l ,v l ) on TU. 

Now, we define T*t m : T*TM -> T*M by (T*T M (a u ),w) = {a u ,w%); 
u,w e T q M, a u £ T*TM and w\ £ T U TM. 

In the following, we recall the construction of the symplectomorphism 
au- To do this, it is necessary to introduce the canonical flip (|12j) on 
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TTM: 



TTM ■ 
TM ■ 



TTM 
-TM, 



as follows: 



where \ '■ 



d 

ds 



d 



Id 



X{s,t) = — 
t=o J as 



d_ 

=odt 



t=o 



X(s,t) 



If 

for 



— > M and \ : M. 2 — > M are related by x(s,t) = x(t,s). 
(g l ) are the local coordinates for M, for TM and (q l , v i , q l , u l ) 

TTM, then the canonical involution can be defined as km uf, v i , q l , v 1 ) = 

In order to describe a a/ is also necessary to define a tangent pairing. 
Given two manifolds M and TV, and a pairing between them (■, •) :MxiV-} 
R, the tangent pairing (•, -) T : TM x TN — > K is determined by 



d 



dt t=o * * at t=o 
where 7 : R -)> M and 5 : R -)■ iV. 

Finally, we can define «m as (aM(^),if} 
TT*M and «; G TTM. In local coordinates: 



t=0 



(7 (<),«(*)) 



(z,km{w)) t , where 2 E 



,Pi,q\Pi) = {q\q\Pi,Pi) ; 

now (q l ,Pi) are the local coordinates for T*M and (q l ,Pi,q l ,Pi) for TT*M. 

The third double vector bundle is T*T*M. The isomorphism : TT*M - 
T*T*M is just given by (3m = ^ui M (see previous subsection). 

Considering the bundles TT*M, T*TM and T*T*M, as well as the sym- 
plectomorphisms olm and Pm, we finally obtain the Tulczyjew's triple 



T*T*M 



13m 



TT*M 



T*TM 



4. Continuous Lagrangian and Hamiltonian mechanics 

In the introduction we have shown how the Tulcyjew's triple is used to 
describe geometrically Lagrangian and Hamiltonian mechanics and its re- 
lationship. In this section we will see that it is also possible to adapt this 
geometric formalism when we introduce constraints into the picture. As we 
have mentioned along the introduction, there are (at least) two methods 
that one might use to derive the equations of motion of systems subjected 
to constraints. We will call them nonholonomic mechanics and constrained 
variational calculus. The classical method to derive equations of motion for 
constrained mechanical systems is the nonholonomic mechanics. The equa- 
tions derived from the nonholonomic methods are not of variational nature, 
but they describe the correct dynamics of a constrained mechanical system. 
In order to obtain the nonholonomic equations, if we have linear or affine con- 
straints, is necessary to apply the Lagrange-D'Alembert's principle. When 
dealing with nonlinear constraints, one should employ the more controver- 
sial Chetaev's rule (see [5j |23] for further details). Since the geometrical 
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implementation of the Chetaev's rule is practically equal to the process in 
the linear case, we shall use it from a pure mathematical perspective. 

On the other hand, the equations of motion of constrained variational 
problems are derivable by using variational techniques (always in the con- 
strained case). These last equations are also known in the literature as 
vakonomic equations. The terminology vakonomic ("mechanics of varia- 
tional axiomatic kind") was coined by V.V. Kozlov ([3],[22j). The main 
applications of the constrained variational calculus appear in problems of 
mathematical nature (like subriemannian geometry) and in optimal control 
theory. 

4.1. Nonholonomic mechanics. A nonholonomic system on a manifold 
M consists of a pair (L, C), where L : TM — > R is the Lagrangian of the 
mechanical system and C is a submanifold of TM with canonical inclusion 
ic '■ C <-} TM. In the following, we will assume, for sake of symplicity, that 
t~m(C) = M. Since the motion of the system is forced to take place on the 
submanifold C, this requires the introduction of some reaction or constraint 
forces into the system. If (p a (q l ,q l ) = 0, 1 < a < n, determine locally the 
submanifold C, then Chetaev's rule implies that the constrained equations 
of the system are: 

d (dh^ dh - \ ^ 



dt \dq* J dq % dq l 

(1) 

< f> a (q i , q i )=0, l<a<n. 

Next, we will describe geometrically the nonholonomic equations. First, we 
need to introduce the vertical endomorphism S which is a (1, l)-tensor field 
on TM defined by 



S : TTM — > TTM 
d 



W, 



dt 



(v x + tTr M (W v J) 



Its local expression is S = <g> dq l . 

If we accept Chetaev-type forces, then we define 

F = S*(TC)°. 

Observe that the vector subbundle F will be generated by the 1-forms fi a = 
|£ dq l because S*(d^ a ) = \x a . 

Now, define the affine subbundle of T^TM given by 

S noh = (dL) o i G + F, 

that is, 

s noh = ^ q i A i^ hh ) eT *TM\ (2) 



dq l dq 1 
<j) a {q,q) =0, 1 < a < n} . 



s 
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Therefore, applying the Tulczyjew's isomorphism olm we obtain the affine 
subbundle 

aj}(Tr h ) = {(q\ Pi ,q\p t )£TT*M\ (3) 
_ <9L 

. _ <9L dcj) a 

Pi — TT~ 7 i A a .,- ) 

oq 1 oq 1 
cf) a (q,q) = 0, 1 < a < n} . 

Define now the nonholonomic Legendre transformation WL noh : C — >• T*M 

by 

FL" o/l = 7r T * M o a" 1 oAoj c . 

The solutions for the dynamics given by (T, noh ^ are curves a : I C R — > 
M such that |(I)cC and the induced curve 7 : R — > T*M, 7 = FL noh (^) 
verifies that ^(/) C a^/ (T 1 noh ). Locally, a must satisfy the system of 
equations ([I]). 

An interesting use of Tulczyjew's triple in order to define Lagrangian 
submanifolds and generalized Legendre transformations within the nonholo- 
nomic framework can be found in |31j . 

4.2. Variational constrained equations. Now, we study the same prob- 
lem but now using purely variational techniques. As above, let consider a 
regular Lagrangian L : TM —¥ R, and a set of nonholonomic constraints 
<fr a (q l , q l ) , 1 < a < n, determining a 2m — n dimensional submanifold 
C C TM. Now we take the extended Lagrangian C = L + \ a (j) a which 
includes the Lagrange multipliers \ a as new extra variables. The equations 
of motion for the constrained variational problem are the Euler-Lagrange 
equations for £, that is: 



d (d(t) a \ dcj) a ~ 
dt V ~W) ~ ~dq T \ ' 

(4) 

cj ) a (q\q i )=Q, l<a<n. 

From a geometrical point of view, these type of variationally constrained 
problems are determined by a pair (C, L) where C is a submanifold of 
TM, with inclusion %q : C TM, and L : C — > R a Lagrangian func- 
tion. Using Theorem 12.11 we deduce that is a Lagrangian submanifold 
of (T*TM,lotm) (see [16J). Now using the Tulczyjew's symplectomorphism 
cxm, we induce a new Lagrangian submanifold cf^ (Sl) of (TT*M, (Itojm)-, 
which completely determines the constrained variational dynamics. Of course, 
the case of unconstrained Lagrangian mechanics is generated taking the 
whole space TM instead of C and an a Lagrangian function over the tan- 
gent bundle L : TM -)■ R. 

Next we shall prove that, indeed, this procedure gives the correct equa- 
tions for the constrained variationalpp dynamics. Take an arbitrary exten- 
sion L : TM — >■ R of L : C — > R, that is, L o %q = L. As above, assume also 



d_ / &L\ dh ■ 90° 
dt \dq l J dq i a dq 1 
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that we have fixed local constraints such that locally determines C by their 
vanishing, i.e: (f) a (q, q) = 0, 1 < a < n where n = 2dim M — dim C. 
Locally 

= {(q\q\^fli)eT*TM\ (5) 
_ <9L d4> a 
^ l dq l a dq l ' 
_ <9L d(f> a 
% dq l ° dq 1 
( p a (q,q) = 0, l<a<n}. 

Observe that locally the conormal bundle i/*(C) = span {d(j) a , 1 < a < n}. 
Therefore, 

a£C£ L ) = €TT*M\ (6) 

_ dh d(j) a 

Pi — TTTJ i A a - —. , 

oq l oq % 
_ <9L d(j> a 

Pi — TT~- i A a - . , 

oq l oq % 
<P a (q,q)=0, l<a<n}. 

The solutions for the dynamics given by a^J (Xl) C TT*M are curves 
7 : I C M -> T*M such that § : I C M -> TT*M verifies that g(J) C 
a^ 1 (Sl)- Locally, if 7(i) = (q l (t),pi(t)) then it must verify the following 
set of differential equations: 

d /asL A ^A_^_ A ^ _ q 

<W,<f) = 0, 

which clearly coincide with equations (HJ). 

Now, we consider adapted coordinates (</*, g a ) to the submanifold C (recall 
that tm{C) = M is now a fibration C — > M), 1 < i < dimM and 1 < a < 
dim M — n, such that 

ic(<z%^) = (?S^,* a (5S?°)). 

This means that (p a (q l ,q l ) = q a — *f/ a (q l ,q a ) = 0. Therefore, we have 

Sl = {(<f,<f,/w*)l (7) 





<9* a 






dL 




^ a ~ dq" 






,<f), l<a<n} 



Observe that (q l ,q a ,fl a ) determines a local system of coordinates for El. 
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Then, 

a£p L ) = {{q\p u q i ,p i )£TT*M\ (8) 



dL d^ 



a 



Pa = T7T- -P, 



Qqd, Qqa 

_ dL d^ a 

Pi — ~ Pa ,• j 

oq l oq % 
q a = V a (q i ,q a ), 1 < a < n} . 

Consequently, the solutions must verify the following system of differential 
equations (see [9]): 

dL d^ a \ dL dV a 



Pa^nr = ir^-Pc 



dt \dq a dq a J dq a dq a 

dL d^> a 



q a = \& a (q\q a ), l<a<n. 
4.2.1. The constrained Legendre transformation. 

Definition 4.1. We define the constrained Legendre transformation FL : 
T,l — > T*M as the mapping FL = tt*m ° 

We will say that the constrained system (L,C) is regular ifFL is a local 
diffeomorphism and hyperregular ifFL is a global diffeomorphism. 

Observe that locally, if as above we consider the constraints q a = *& a (q l , q a ) 
determining C, then 

dL d^ a 
FL(q\q a ,jx a ) = (q l ,p a = — - fi a -^-,p a = fi a ) ■ 

The constrained system (L, C) is regular if and only if (^§^ - /"a^prf^) 
is a nondegenerate matrix. 

Next, define the energy function El '■ £i — > K by 

Ehiau) = (at u ,v%) - L(u), a u eT, L ,u £ C = z'c(C) 

Locally, we have 

dT d^> ct (n i n a ) 
E L (q\ q a , M = q a Q~ a ~ K f^- 1 f + ~^ a {q\ <f) ~ L{q\ «f ). 

Remark 4.2. The constrained Legendre transformation allows us to de- 
velop a Lagrangian formalism on E^. Indeed, we can define the 2-form 
uil = (WL)*ujm on Si and it is easy to show that the equations of motion 
of the constrained system are now intrinsically rewritten as 

ix^L = dE L . 

In consequence, we could develop an intrinsic formalism on the Lagrangian 
side, that is a Klein formalism (|12|. |17j . |21| . |25j ) for constrained systems 
without using (at least initially) Lagrangian multipliers. Moreover, notice 
that the constrained system is regular if and only if ljl is a symplectic 2-form 
on S L 
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Then, if the constrained system (L, C) is hyperregular we can define the 
Hamiltonian function H : T*M — > R by 

H = E L o (FLy 1 , 

and the corresponding Hamiltonian vector field Xh by ix H ^M = d.H. In 
this particular case we have that 

lmX H = X H (T*M) = p-\dH(T*M)) = a" 1 (S x ) . 

The second equivalence, Xh(T*M) = cf^ will be studied in § 14.41 

The next diagram summarizes the above discussion: 



TT*M ^ T*TM 




In the singular case, it is necessary to apply the integrability algorithm to 
find, if it exists, a subset where there are consistent solutions of the dynamics 
(see [S], pj). 

4.3. Comparison of nonholonomic and variational constrained equa- 
tions. Continuous picture. Let consider a system defined by the La- 
grangian function L : TM — > M. and an independent set of constraints 
(f> a (q l ,q l ) = determining the submanifold C C TM. 

As shown in § 14.11 and § 14.21 the solutions of the nonholonomic dynam- 
ics are geometrically described by the submanifold Y 1 noh =— )■ TT*M, while 
the solutions of the constrained variational dynamics are given by the La- 
grangian submanifold > TT*M, where L = L| c : C — > R. 

Our aim is to know whether, given a solution of the nonholonomic prob- 
lem, it is also a solution of the constrained variational problem. In order 
to capture the common solutions to both problems, we have developed the 
following geometric algorithm. Consider the fibered product T*M © T*M, 
where we choose the local coordinates (q l ,pi,7Ti); consider also the tangent 
bundle T(T*M © T*M) which can be identified with TT*M TT*M, 
which fibers over TM. Under these considerations, construct the submani- 
fold S cons «-»■ TT*M © T TT*M as follows: 

T7T M 

S cons = { [X at ,Y Pt ) G TT*M® t ^ m TT*M/T aq ir M (X aq ) = T^ M (Y 0q ), 
for ^er'j^eEi}. (9) 

It is quite clear that the submanifold £ cons gathers together both nonholo- 
nomic and constrained variational dynamics. Locally, Y, cons is determined 
by the coordinates (q l ,pi,iTi,q' l ,pi,7ri) obeying the nonholonomic and con- 
strained variational conditions respectively presented in ([2]) and ©, that 



12 



M. DE LEON, F. JIMENEZ, AND D. MARTIN DE DIEGO 



IS 



Pi 



Pi 



8L 
dqi- 



— Qqi 



c)q' 



dg 



- 4_ X ot P 

a <9q* ' 



7T; 



dh , d<j><* 



subject to 4> a {q l , q l ) = 0. Here, X a and fx a are the nonholonomic and vari- 
ational constrained Lagrange multipliers, respectively. Finally, in order to 
find the common solutions we will need to apply the integrability algorithm 
described in [9] [32] . 

The following diagram shows the bundle relations: 



TT*M < 



TT*M ■ 



TM 



( T T* M < T T* m) 



T*M T*M ■ 



M 



KM 



where Ttt m : TT*M e T7r TT*M — > TM and -\, : T*M \ T*M — > M 
denote the fibrations of the Whitney sums over TM and M, respectively. 

As a simple example, consider the case of linear constraints, namely 
4> a (q i ,q l ) = q a — f2(q l )q a = 0. The relationship determining £ cons pre- 
sented in ([25]), as well as the integrability algorithm, implies the following 
equation: 



where 



R a 



9<P° b 



ah 



q y Q R% b = o, 



(10) 



dq a dq b ' ra dqP 

can be considered as the curvature of the connection T in the local projec- 
tion p(q a ,q a ) = (q a ) such that the horizontal distribution H is given by 
prescribing its annihilator to be 



{dq a - ip«dq a , 1 < a < m}. 



See more details in [S]. 



4.4. Lagrangian and Hamiltonian mechanics relationship. In this 
section, we shall discuss the converse case, i.e., starting from a Hamilton- 
ian system we shall show that it is possible to construct a constrained La- 
grangian system providing the same dynamics. 

Since ttm '■ T*M — > M is a vector bundle, then it is possible to define 
the dilation vector field or Lioville vector field A* G £(T*M), which is the 
generator of the one-parameter group of dilations along the vertical fibres 
fiq — > e t fi q , fj, q G T '* Q . The dilation vector field is locally expressed by 
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Given a Hamiltonian function H : T*M — > R, the Hamilton's equations 
are written in canonical coordinates by 

■i 9H 

(11) 

• - dH 
^ % dq l ' 

The solutions of the Hamilton's equations are just the integral curves of the 
Hamiltonian vector field given by 

ix H u M = dH, (12) 

where ujm is a symplectic form on T*M . 

Given the Hamiltonian function, one defines a function FH : T*M — > 
TM, i.e. the fiber derivative of H (see [I]), by 

{WH(a g ),0 g ) = — H(a g + t[3 q ), 



dt 



t=o 



where both a q ,(3 q 6 T*M. In local coordinates, 

f) H 

VH{q\ Pi ) = {q\—). 

OPi 

Assume that the image of T*M under FH defines a submanifold C of 
TM. Mimicking the Gotay and Nester's definition ([14|.[15]). we implicitly 
define the function L : C — > K by 

L o FH = A*H - H. (13) 

The function L : C — > M. will be well-defined if and only if, for any two points 
a q , p g G T*M such that FH(a q ) = FH(f3 q ), we have that (A*H - H)(a q ) = 
(A*H — H)((3 q ). Obviously, without additional assumptions there is no 
reason why this should be true. The following definition states under what 
conditions such projection L exists. 

Definition 4.3. A Hamiltonian H : T*M — >■ R is almost-regular if 

FH(T*M) = C is a submanifold of TM and FH : T*M -> C C TM is 
a submersion with connected fibers. 



T*M — ^ TM 

FH(T*M) = C 

Under the assumption that the Hamiltonian H is almost-regular, it is only 
necessary to show that expression (|13p defines a single-valued function L : 
C — > R, or, in other words, that A*H — H is a constant function in the 
fibers. Since each fiber of the submersion FH is connected, it is sufficient to 
consider the infinitesimal condition, i.e. to show that 

C Z {A*H-H) = 0, for all Z E ker(F^) , 
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where Lz is the Lie derivate in the Z direction. Working in local coordinates, 
Z will be of the form 

d d 2 H 
Z = Zi—, with Zi - — — = 0, for all 1 < j < n . 

dpi dpi dpj 

Since FH*(Z) = 0, the last condition can be obtained taking into account 
that (a, FH*(Z)) = (FH*(a) , Z) = 0, for a an arbitrary point of T*TM 
such that ktm(o-) = ttm (FH*(Z)). Then 

dpi 

dH d 2 H dH 

dpi dpidpj dpi 

In what follows we will assume that H satisfies the almost regularity prop- 
erty. 

Theorem 4.4. The following equality holds 

a M {X H (T*M)) = S L , 
where ocyi is the Tulczyjew's isomorphism. 

Proof. Take W x G Z Xh = X H (T*M) C TT*M and take a M (W 1 ) G T*TM. 
We need to prove that 

(a M (Wi),f/) = (dL,U), 

for all U G TC C TTM, such that ttm(U) = Titm(Wi). This is equivalent 
to the equality 

(a M (W 1 ) , WH*(W 2 )) = (dL , FH*(W 2 )) , (14) 

for all W 2 G TT*M such that t T m (Fi/*(W 2 )) = Ttt m {Wi). Therefore, 
regarding (TT3]) the previous equality turns out to be 

(FiJ)* (a M (Wi)) = (FH)* (dL) = d(A*H - H) , (15) 

where the right hand of the equation comes directly from (|13p . If locally 
Wi = q 1 ^ +Pi^-, in other words W\ = (q l , pf, <f , p\) , then 

(FHna Mm) = {^V^)dq^p^ dPl 

d 2 H dH\ { | d 2 H d 

^ dpjdq 1 dq % ) ^ ^ J dpjdpi l ' 

where we consider pi = — since we are dealing with a Hamiltonian vector 
field Xh (see equations (fTTj) ). From condition, "Prj\/F£/*(W2) = ^m^i) 
we also deduce that q 1 = |^ . But this is true since W\ G T,x H ■ 

On the other hand, a straightforward computation leads to check that 
d(A*H - H) is exactly flUD. □ 
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5. Discrete equivalence 



5.1. Discrete nonholonomic mechanics. A discrete nonholonomic sys- 
tem is determined by three ingredients: a discrete lagrangian : M x M — > 
R, a constraint distribution PonM and a discrete constraint submanifold 
Cd C M x M with canonical inclusion ic d : Cd <->• M x M '. Notice that 
discrete mechanics can also be seen within this case, where T> = TM and 
Cd = M x M. Notice also that, in the discrete version of Lagrangian me- 
chanics, the tangent manifold TM is substituted by the cartesian product 



Define the affine submanifold cT*(Mx M) by Y%° = (dL d ) oi Cd + 
Fd, where is the vector subbundle of Tq (M x M) given by 



Here, pr\ : M x M — > M is the first projection onto M. 

The symplectic manifold (T*(M x M) , wjuxm) is symplectomorphic to 
(T*M x T*M, ft), where ft = prJ^M - pr^ M and pr i : T*M x T*M -> 
T*M, i = 0, 1, are the natural projections of T*M x T*M onto T*M. The 
symplectomorphism is given by 



where {qo,qi) £ M x M. 

Using T we induce the affine subbundle T(S^°' 1 ) of the symplectic mani- 
fold (T*M x T*M, ft). The dynamics is then determined by the sequences 



750,7,1, • ■ ■ ,7«r ™ch that (7«,7« +1 ) e T(S^), < i < A - 1 (see [2D]). 



We will describe now the equations in terms of local coordinates. As- 
sume that Cd is defined by the vanishing of the following set of independent 
constraints: (j>d(lo,Ql) = 0, 1 < a < n, where n = 2dim M — dim C and 
V° = span{w a = uif dq 1 }. Therefore, 



M xM. 




T:T*(MxM) T*M x T*M 

7( ?0 ,?i) = (7 go >7gi) ^ (-7go>7?i) 



noh 



(Mo)i,(m)i) er*(MxM) 

(A*o)i = ^"T + ( A o) a ^i , 



(17) 




0d(?o,9i) = O, l<a<n}, 



where (A Q )o are Lagrange multipliers to be determined. Then, 



fx^noh\ 
K^d ) 



M,(po)i,ql(Pl)i)£T*MxT*M 



(18) 




0d(?o,9i) = O, l<a<n}. 
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The solutions of the dynamics are therefore given by the equations 

dq{ dq\ 
<t>d{qk,qk+i) = 0. 

These equations are traditionally written in the following manner 

0d{Qk,qk+i) = 0, 

which are the expression of the discrete nonholonomic equations (see [TO] 
for more details). 

5.2. Constrained discrete Lagrangian mechanics. A discrete constrai- 
ned system [29] is determined by a pair (Cd, Ld) where Cd is a submanifold 
of M x M, with inclusion %c d ■ Cd ^ M X M, and Ld : Cd — > K a discrete 
Lagrangian function. 

Using again Theorem l2.1l we deduce that is a Lagrangian submanifold 
of (T*(MxM), umxm)- 

Using T we induce the Lagrangian submanifold T(S^ d ) of the symplectic 
manifold (T*M x T*M,9) (see the diagram below). 

T*(M x M) ^ T*M x T*M 



^L d -T(S id ) 

The dynamics is determined by the sequences j qo , r y qi , . . . , 7 9JV such that 
(7 9i ,7g i+ i) G T(S L J, < i < N - 1 (see [201 [29]). Observe that 

7 9fe G T 9 * fc M npro(T(S L J) n m (T(S L J), 1 < A; < JV - 1. (19) 

After determining intrinsically the dynamics, as in the continuous case, we 
now consider local expressions. Take an arbitrary extension L^MxM-fR 
of Ld : Cd — > M, that is, o = L^. Assume also that we have fixed 
local constraints such that determines the submanifold Cd- This definition is 
performed by the vanishing of the following set of independent constraints: 
4>d(.l0i 9i) = 0, 1 < a < n where n = 2dim M — dim C. 

Locally 

Si, = {(AiWi,Wl)6T*(MxM)| (20) 
( w)i = fS + (A l)a 8 « 



^d(9o,9i) = 0, 1 < a < n} , 
where (Ai) Q are Lagrange multipliers to be determined. 
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Therefore, 

T(£ L „) = {(4>(Po)i,<A,(pi)i)eT*MxT*M\ (21) 

/_ x _ dh d V 

dql dql 

<t>d(lo,qi) = Q, l<a<n}. 

The solutions of the dynamics come from (|19p and are given by the equations 

9hd d<j>% dhd d(f>°l . 

-rr-r(qk-i,qk) + {^k)a^-f{qk-i,qk) = -jr-f{qk,qk+i) - (Afc + ij a —- j-{qk,Qk+i)) 
aq 1 Oq 1 oq Oq Q 

<i>d(qk~i,qk) = o, 
0d((?fci?k+i) = o. 
These equations are traditionally written as 
D2^d{qk-i,Qk) + Dih d (q k ,q k+1 ) 

+{\k) a D 2 4>2(qk~i,qk+i) + (h+daDi^iqk^qk+i) = o, 

which are the expression of the discrete vakonomic equations (see |5j for 
more details). 

Now, as a particular case, we assume that we can choose adapted coordi- 
nates (qQ,qf), 1 < i < dimM and 1 < a < dimM — n, on C<2 in such a way 
the inclusion is written as 

icM,ti) = (<&,&*d(4,ti))- (22) 

In other words, we can write the constraints as <fi°[(qo, qi) = qf — ^(^oj 9?) = 
0. Thus, locally we have that 

Sl, = {(<Zo> Qii (jM))i} (Mi)i) e T*(M x M)| (23) 

tf = Wo.9i). 1<«<^}- 

Observe that (q^qf, (Mi)a) gives a local coordinate system for 

Therefore, we obtain the following expression of the Lagrangian subman- 
ifoldT(S L J: 

T(S L J = {(<t ,fa>)i,qi,(pi)i)€T*MxT*M\ (24) 



0$ dgf 



o 
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Consequently, the solutions must verify the following system of difference 
equations: 

(D 2 ) a (L d - (pjfc)o*3)(?fc_i,?fc) + {Dx) a {L d - (Pfc+i)«*d)(gfc,gjfc +1 ) = 0, 
(Pk)/3 + (Di)pL d (q k ,q k+1 ) - {Pk+i)a(D 1 )p^(q k ,q k+1 ) = 0, 

-l), 

where (Dj)i just means being j = {1, 2} and / = {a, a}. 



5.2.1. The constrained discrete Legendre transformations. 

Definition 5.1. [29J We define the constrained discrete Legendre transfor- 
mations FL^ : Y,L d — > T*M as the mappings 

WL- = pr o(T)\ , 



FL+ = prio(T) 



We will say that the constrained system (L d ,C d ) is regular if¥L d is a local 
diffeomorphism and hyperregular if¥L^ is a diffeomorphism. 

Remark 5.2. It is easy to prove that FL^~ is a local diffeomorphism if and 
only if FL^~ is a local diffeomorphism; therefore, it is possible to character- 
ize the regularity of the constrained system using any of the two Legendre 
transformations (see |29|). 

Observe that if we consider the local constraints qf = ^^ilb^Qi) deter- 
mining Cd, then 

FLj(g*,(tf,0*i) a ) = ( g * J (p ) i = --^ + ( Ml ) a _f) 

°% °% 

= -(£>i)<(£ d - (A*i)«*5)(flb,ffi)) • 

FL+tf ,fiM a ) = (q a 1 ,q? = * a M, q a 1 ), 

. , dL d ... 

ipi)a = (D 2 )a(L d - /i«*d)(gb, gi) , (Pl) a = Ol)a) • 

So, the constrained system (L d , C d ) is regular if and only if the matrix 
(A- A- ) - (^L _ u i dK 



^dq l dqf Q dq^dq^ dq l Q 
is nondegenerate. 

Remark 5.3. The constrained Legendre transformations allow us to define 
a univocally presymplectic 2-form on 

u Ld = (WL-yujM = (¥L+)*u M 

which is symplectic if the system is regular. 
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Then, if the constrained system (Ld, Cd) is hyperregular, we can define 
the discrete dynamics determined by T _1 (Si d ) as the graph of the canonical 
transformation (FL+) o (¥L d )~ l : T*M -> T*M, that is, 

Graph (FL+)o(FL,)- 1 = T(S L J. 

5.3. Comparison of nonholonomic and variational constrained equa- 
tions. Discrete picture. Let consider the system defined by the discrete 
Lagrangian function : M x M —> R and a set of constraints ^(qo, qi) = 
determining the submanifold Cd C M x M. 

As shown in § 15.11 and § 15.21 the solutions of the discrete nonholonomic 
dynamics are geometrically described by the affine subbundle T(S^°' 1 ) C 
T*M x T*M, while the solutions of the discrete constrained variational 
dynamics are given by the Lagrangian submanifold Y(Y,L d ) C T*M x T*M, 

where Ld = LJ „ : Cd — > R- 

d 

Given a solution of the discrete nonholonomic problem, we want to know 
when is also a solution of the associated discrete constrained variational 
problem. To capture the set of common solutions to both problems, we 
have developed the following geometric integrability algorithm similar to 
the continuous one. 

As in the continuous case, take the Whitney sum T*M © T*M, and the 
cartesian product 

(T*M © T*M) x (T*M © T*M) = (T*M x T*M) ®n M xn M (T*M x T*M). 

Construct the submanifold £™ ns {T*M © T*M) x (T*M © T*M) as 
follows: 

Sr s = { (7 90 , 7 91 , 7 90 , 7«n) 6 (^*M x T* M) TMXffM (T*M x T*M) / 

(Tft.T^esf.MeSU- (25) 

It is quite clear that the submanifold $]^ ons gathers together both non- 
holonomic and constrained variational dynamics and applying the discrete 
version of the integrability algorithm developed in |20| we will find the set 
where there are common solutions for both dynamics. 
The following diagram shows the bundle relations. 

>:;/"'- • (T*M © T*M) x (T*M © T*M) M x M 

■km x km 



T*M © T*M — + M 

^M 

where the vertical arrows represent the projections onto the first and second 
factor of the respective cartesian products. 

5.4. Construction of a discrete constrained problem from a Lagran- 
gian submanifold of T*M xT*M. Given a Lagrangian submanifold A of 
(T*M x T*M, $7) we will construct, under some regularity conditions, a 
constrained Lagrangian problem given by a submanifold Cd C M x M and 
a function : Cd — > R. 
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We first construct the Lagrangian submanifold Y _1 (A) of (T*(M x M), 
^Mxm) and we assume that the restriction of ujmxM to Y _1 (A) is exact, 
that is, we have a generating function S : T _1 (A) — > R. In addition, we 
suppose that the image of Y -1 (A) by ttmxM is a submanifold Cd of M x M, 
and, also, that (7Tmxm)|y-i(a) ^ s a submersion with connected fibers (see 
the diagrams below). 

T*(M x M) > T*M x T*M 

T-HAH ^1 A 



T*(M x M) 




Theorem 5.4. Under the previous conditions the function S : T X (A) — >■ M 
is (7TMxM)| T -i( A )-proiecto6Ze onto a function La : Cd —t K. Moreover, the 
following equation holds 

T- 1 (A) = S Ld . 
Proof. The submanifold T _1 (A) is defined as 

T- X (A) = {( 7go ,7 9l ) G T* {qom) {MxM) j i* r _ 1(A) 0MxM(l qo ,l qi ) = dS( 7qo , lqi )}. 

(26) 

By definition of the Liouville 1-form OmxM we have that 

(OmxM , ker T itmxm) = 0. 
By applying the chain rule T(ir M xM ° S-i (A) ) = Ttt M xM ° Ti T -i (A y is 
easy to check that T i , (kerT7TMxM| I C kerTyr^xM- Thus, we 

(A) y 'T-i(A) y 

finally deduce that S is projectable into Ld : Cd~> K, i.e. 

5" = (7tmxm| (27) 

't-!(A) 

since 7Tmxm(T _1 ( a )) = C<i- 

On the other hand, we can define a tangent vector X = {X lqo , X Jqi ) G 
TT _1 (A) C TT*{M x M), such that T^^ttmxmPO = G 

T (qo,qi) C d, b y 

( (7 9o . 7«i ) , Ko . )) = ( d5 (7<?<> . 7«i ) , -X") , 

where (7 ?0 ,7 ?1 ) G T _1 (A). The equation above comes directly from the 
definitions of both the Liouville one-form and the Lagrangian submanifold 
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T X (A) in (f26jh Regarding equation ([27j) and taking into account that the 
pullback and the exterior derivative commute, we arrive to 

(dSKoiTgi) ' x ) = ( dL d(qo,qi) , T7I mxa-/| t _ 1(a) PO) 

= (dL d (q ,qi) , (u qo ,u qi )). 

In the last line of the expression just above, we recognize the definition given 
in Theorem 12,11 of , that is 

^L d = {(lg ,l qi ) G T ( qo , qi )( M X M ) I ((7<?o>7 gi ) , (Ug ,U qi )) = 

= (dL d (q ,qi) , (u qo ,u qi )) for all (w 3o ,tt gi ) G ?( 30)9l )C d }' 
In consequence, we deduce that T _1 (A) = S Ld . 

□ 

Remark 5.5. We would like to point out that the proof of Theorem 15, 41 can 
be easily extended to the continuous case. Namely, let consider A C TT*M 
a Lagrangian submanifold. Under some regularity conditions, it is possi- 
ble to construct a constrained Lagrangian problem given by a submani- 
fold C C TM and a Lagrangian function L : C — > R. If we consider the 
Lagrangian submanifold ajvf(A) C (T*TM , ujtm), where ajvf is the Tul- 
czyjew's isomorphism, we can build an analogy with the discrete case by 
assuming that the restriction of ojtm to «m(A) is exact, that is, we have a 
generating function S : a a/ (A) — > R; 

«Af (A) = { 7 G T^TM I £ m(A) 0tm (7) = d£( 7 )} , 
where v q G TM such that T M (v q ) = q G M. In addition, we suppose 
that 7T TM ^aj\/(A)^ is a submanifold C C TM, and that 7I "ta/| Qm (a) i s a 
submersion with connected fibers. 

Again, by the definition of the Liouville one- form Qtm w e have that 
(0 TM , kerT7r TM ) = 0, and consequently that S is projectable into L, that is 

S = ( 7I "tmL m (A))* 1 '- 

Following similar arguments that in the proof of Theorem 15.41 we deduce 

that a a/ (A) = S^. 

Since we are not fixing the Lagrangian submanifold A, this is a more 
general result than that one provided by Theorem l4.4l i.e. aM(Xjj(T*M)) = 
Til. Nevertheless, among all the Lagrangian submanifolds A of T*TM, we 
choose Xh(T*M), that is the image of the cotangent bundle T*M by the 
Hamiltonian vector field provided by the equations ix H ^ = dH, in order 
to stress the relationship between Hamiltonian and constrained Lagrangian 
systems. 

6. Examples 

6.1. Linear Constraints. Let consider a dynamical system denoted by the 
Lagrangian 

Hv q ) = \ 9(v q , v q ) - V(q) = ^(f (f - V(q) (28) 
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where v q € T q M, with local coordinates v q = (q l ,q l ), g is a Riemannian 
metric with components (gij). Moreover, V : M — > R is a potential function. 
Additionally, the system is subject to the linear constraints 

<p a (v q ) = q a -K(q)q a , (29) 

where q % = {q a , q a }. Locally, the constraints define a submanifold C C TM. 
Moreover, we have the restriction of L to C, L : C — > R. In local coordinates, 

L{q\q a ) = \l ab q a q h -V{q), 

where 

lab (q) = g ab + gac£t(q) + 9baK{q) + g a pK{q)^ b {q)- 

Observe that (jab) is invertible since g is a Riemannian metric. 
Using expression (|S]), we can find local coordinates for (£l): 

Pa = lab{q) q b ~ Pa^aiq) > 

" " I-Wft-Bf-*??®*- 
q" = r«(?)««. 

The Legendre transformation is defined by FL = tt* m (a^ 1 ) L , or locally 
by: 

F£(g*, g a , Aa) = (9*, 7a&g° - Aa 1%), Aa) • 

Since f g^s - Aa^^t) = (7a&) the constrained system (L, C) is regular. 
Moreover, the energy function ^ : El — > R is precisely 

W,g a ,Aa) = ^7a6g a g 6 -ng)- 

Therefore, the Hamiltonian function can be expressed by H = El o (FL) -1 : 
T*M -)• R 

^,p) = ^7 a6 (g) J Pan + U(«), (30) 
where 7 a6 7 6c = 5 c a and P a = p a +p a T a t (q). 

6.1.1. Discretization: symplectic Euler method. Taking into account equa- 
tions ([28]) and (f29|) , we define the discrete Lagrangian L^MxM-jB and 
the set of independent constraints in the following way (see [30J for more 
details): 

M<7o,<?i) = hh(q , q% q ° ) = ^rgij(qo){q\ -<?o)(<?i - Qo) -hV(qo), 

(31) 
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Next, we can explicitly obtain the coordinates for the submanifold T(S^ d ) 
given in equations (|2ip . namely 

(Po) a = j^9aj{qo)(q{ ~ 4) + hd °y (lo) - ^da9ij(qo){q\ - ?o)(«i - 4) 

+ (Ai)^ d a T^q ) {q\ - ql) - (X^ a T a a (q ) , (32) 

(Po) a = ^5«i(9o)(9i - 4) + hd <*V (go) - ^d a g lj (q Q ){q l 1 - q l )(q{ - <? Q ) 
+(Ai)„ + (Ai) / ,0 a r£(< to )(gJ-gg), (33) 

(pi), = ^ai(gd)(«i-^)-(Ai)ars(«,), (34) 

(Pl) a = J l 9aj(qo)(4 - 9o) + 

(??-«?) = C(g )(g?~ go a ), 

where 3 a , (9 Q mean and respectively. It is important to note 
that (|32|) is a set of 2m + n equations with 2m + n unknowns, which are 

(qi) a ,(qi) a ,(p0o» (pi)o and ( Al )«- 

Alternatively, we can apply the so-called Euler symplectic method (see 
Pi=Po-h—{q Q ,Pi), qi = qo + h—{q Q ,pi) (35) 



to the Hamiltonian function H : T*M — t M defined in (|30p . We deduce the 
following set of equations: 

(9i)° - (gor + ^^i),, 

( gi ) a = («,)" + ftrjy 6 ^, 

(36) 



(pi), = (po) a -^ a y-^|i(9 Q7 fcc )(F 1 ) b (p 1 ) c + 7 bc (p 1 ) h a a (p 1 ) ( 
(pi)a = M Q -^a^-^{i(a Q7 ab )(A) a (Pi) h + 7 Qb (A) Q 9 Q (Pi) 

where (P x ) a = {pi) a + (pi) a r°, and V, 7 afe , r„, (P\) a are evaluated at go- 
Regarding equations (|32p . is easy to express Ai in terms of p and q. Since 



= cWS = d a (7 cb 7 M J = ^a7 c "J 76d + T (d a %d) , (37) 

and (9 a 7 c6 ) 7bd = —J cb {dajbd) , is easy to check, after a straightforward 
calculation, that equations ([32]) reduce to (f36j) . 



6.1.2. Discretization: Midpoint rule. Define the discrete Lagrangian and the 
discrete constraints using the midpoint rule (see [30] for more details), that 
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IS: 



Mgo, gi) = WL( — - — , — - — ) = 

^(^y^ )(<?1 - qoY(<n qo) j hv( 



qi + qo, 



(38) 



gl ~ gQ A pq/ gl + g0 W gl - g0 

Now, we can explicitly obtain the corresponding equations derived from 
the submanifold T(Ex d ). A straightforward computation shows that these 
equations are equivalent to the corresponding ones derived from the so-called 
midpoint rule 

dH f q 1 + q Q P!+p \ dH (q x +q Pi + p \ , , 

* = * + V (—'—J' Pl=P °-^(^-'^— J' (39) 
which is a symplectic method of order 2. 



Vrf 



T*M 



•pd 



T*M 



(40) 



£ L - - ^ E L , 



(4>H)d 



T*M < T*M 
T*M^--T*M 



(<f>H)d 



(41) 



T*M 



T*M 



•pd 



The results of the previous examples can be summarized in the diagrams 
(|40p and (|4ip . which are explained in the following lines. 

• Prom right to left: (pn is the Hamiltonian flow derived from (|lip 
applied to a Hamiltonian function H, (pd is the discretization of that 
flow (concretely the symplectic Euler method (|35p and the midpoint 
rule ([39]) ). {4>n)d is the discrete flow in T*M provided by (pd- 

• From left to right: (f>L is the constrained Lagrangian flow resulting 
from equations ([U applied to the continuous Lagrangian, ipd is the 
discretization applied to that Lagrangian ( ([31"]) and ([38])), 4>L d is the 
discrete flow within l[23|) due to ipd, (j>j L ± is the discrete flow in 
T*M x T*M dMD- 

In order to be more explicit, ipd '■ T*M — >• T*M represents the discrete flow 
generated by applying a symplectic method to the Hamiltonian equations 
(e.g. ([35]) and ([39]) ). On the other hand, ipd represents the discretization 



mapping, that is (pd ■ TM 

Mq,q) = ( qj ^, 9 - h± jF M )- 



M x M, e.g. <pd(q,q) = (qk 



or 
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As expected from Theorems 14.41 and 15.41 what we explicitly show is that 
4> ¥L ± = (4>H)d using some discretizations (we have depicted the particu- 

d 

lar cases of the symplectic Euler methods and the midpoint rule). In other 
words, the diagrams (|40p and (|4ip are commutative in those particular cases, 
as also is when the discretization 0^ corresponds to the exact discrete La- 
grangian's (see [30]). Therefore, using a discrete variational integrator for 
the constrained continuous Lagrangian system or applying a symplectic in- 
tegrator to the associated continuous Hamiltonian problem are equivalent 
approaches. 

6.2. The Martinet case: symplectic integrators for sub-Riemannian 
geometry. Let us consider the Hamiltonian function H : T*R 3 — > R 

*<».ri- 5 prf^) 1 (42) 

where q = (x,y,z) T £ R 3 and p = (p x ,P y ,Pz) 6 (R 3 )* — R 3 - From (|4~2j) we 
can locally define Xh(T*M), particulary Xh(T*R 3 ), through the Hamilton- 
ian equations, i.e: 

• - _l vL ■ - Pvl 

X — Px-\-Pz 2 , Px — (I+/3~rp> 

y = (l+B x)* ' Py = ~ (px+Pz^) p z y, (43) 



(l+/9a;) 

Z = [Px+PzK- ) Pz = 0. 



The associated Legendre transform FH is in this particular case written as 



v 2 , p y ( , y 2 \ y 2 



FH(x,y,z;p x ,p y ,p z ) = (x,y,z; (jp x +p z —), / 1 + p x y > (P*+P*y 1 y)- (44) 
Looking at (|43p and (|44p is easy to realize that 

C C TR 3 = / (x, y, z ; x, y, z) s.t. z = y- x 

Next, we will obtain the Lagrangian function L : C — > R using the implicit 
equation given in (fl~3j) : 

,2~ 




Finally, using FH we arrive to 

L(x, y, z, x,y) = - {x 2 + (1 + (3 x) 2 y 2 ) . (45) 

Consequently, our approach allows us to conclude that the Hamiltonian sys- 
tem (|42p is equivalent to the Lagrangian one (|45p subject to the constraints 

2 

z = \ x. We clearly recognize in (|45p a Martinet sub-Riemannian structure 
(0, [S], [Z])) which is described by the triple (U , A , g). In this triple, C7 is 
an open neighborhood of the origin in M 3 , A is a distribution corresponding 

2 

to A = ker a for a = dz — H^-dx and g is a Riemannian metric. In the 
particular case g = dx 2 + (1 + /3x) 2 dy 2 , we deduce that (|4"5"|) corresponds 
to L = h\ c where h(q, q) = \ g , ^), where ^ = +y^ + z-j%. In 

addition, the constraints are given by q ( J^) = 0. 
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6.2.1. Discrete Case. Let us consider the symplectic Euler method (|35p . Is 
easy to see that a type-2 generating function of the approximated Hamil- 
tonian flow is 

H + (q ,Pi) = qoPi + hH(q ,p 1 ), (46) 
where H is the Hamiltonian function H : T*M 3 — > M. In other words: 

dH+(q , Pl ) dH+(q , Pl ) 
oqo dpi 

Under these considerations, we can define the local coordinates for Y -1 (A): 

T (A) = -U , , ~ , Pi 

{ dq dp x 

Now, projecting T _1 (A) onto M 3 x M 3 , we obtain that 

( y 2 

x\ = x + h Upi) x + (pi) z -^- 

Vl = yo + h --^ , (47) 
(1 + /3x ) 2 

zi = *o + /*(W + (pi)*f)f ■ 

v 2 

From the last equations we obtain the constraint {z\ — zq) = -£-(x\ — xq), 
which defines the submanifold 

y 2 

C d = {(a: 0) yo,«o;»i 3 yi,«i) G M 3 x M 3 | (z 1 - z ) = -j( x i - x )}. 

The next step to completely determine the discrete constrained Lagrangian 
system is to obtain Ld- In that sense, we take the usual generating func- 
tion S mentioned in Theorem 15.41 Transforming S into a type-2 generating 
function (|24J) we arrive to the implicitly defined expression: S(qo,Pi) = 
pi qi — H + (qo,pi), which, taking into account (|4"o]) leads to 

S(q ,px) = h (pi dH{ !£ , i Pl) ~ H(%,Pi)) ■ (48) 

We have shown in Theorem 15.41 that S is ttmxm\ t -i^ -project able onto 
Ld '■ Cd — » M. Thus, from (|48p and according to equations (|4"7|) we finally 
arrive to 

,2 



L d (q ,qi) = ^ ((xi - x ) 2 + ^-(yi - y ) 2 ) . ( 19) 



for 



This is the expected result as is easily seen taking into account the contin- 
uous Lagrangian (|4"5j) and the submanifold C C TM defined by the contin- 

2 

uous constraint z = ^- x. If we define both the discrete Lagrangian and the 
discrete submanifold Cd in the usual symplectic Euler discretization Q30J), 
i.e. Ld(qo,qi) = hL(qQ, qi ^ qo ), we easily obtain P9|) and (|5T)|) . 
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6.3. Symplectic Stormer-Verlet method. Due to its well-behaved fea- 
tures, namely reversibility, symplecticity, volume preservation and conserva- 
tion of first integrals, the Stormer-Verlet method is one of the most impor- 
tant examples in geometric numerical integration (see [18] and references 
therein). For a Hamiltonian system determined by H : T*M — > R, the 
Stormer-Verlet method reads 

Pk+i/2 = Pk - ^H q (q k ,p k+1 / 2 ), 
Qk+i = Qk + 7; {H p (q k ,Pk+i/2) + H p (q k+1 ,p k+1/2 )) , (51) 

Pk+i = Pk+i/2 - -^H q {q k+ i,p k+l i 2 ), 

where q k € R n , p k £ (W 1 )* and H q , H p are the derivatives of the Hamiltonian 
function respect q and p, respectively. As we did for the momenta, we can 
fix an intermediate configuration point q k+ \/ 2 = qk + ^H p (q k ,p k+ i/ 2 ) and 
consider (f5"Tj) as a two step integrator: 

h . . 
Qk+1/2 — Qk + -^^p{Qk,Pk+i/2), 

(52) 

Pk+l/2 = Pk~ -^H q (q k ,p k+1 / 2 ), 



Qk+i — Qk+i/2 + -^H p {q k+ i,p k+ i/ 2 ), 

(53) 

Pk+l = Pk+l/2 - ^lilk+l, Pk+l/2)- 

Equations (|52p and (|53p show the well-known fact that the Stormer-Verlet 
method is the composition of two different symplectic Euler schemes. In 
addition, it is clear that they are respectively generated by the 2- and 3- 
type generating functions 

H + (q k ,p k+1 / 2 ) = p k +i/2Qk + ^H(q k ,p k+1/2 ), 

H~(q k +i,p k+ i/ 2 ) = P k +i/2lk+i - -^H(q k+1 ,p k+1 / 2 ), 
which, taking into account that p\ &q\ — po dqo = dS(qo, q\), lead to 

S + (q k ,Pk+i/2) = Pk+i/2Qk+l/2 ~ Pk+l/2 ~ -^H(q k ,p k+1 / 2 ), (54) 
S~{p k +l/2,qk+l) = -Pk+l/2 fffe+l/2 + P/c+l/2<Zfc+l - ^H(q k +i, Pk+l/2)- (55) 

Now, as shown in [23], we can construct a 1-type generating function 
S{qk,Qk+i) by 

S{qk,Qk+i) = S + (q k ,p k+1/2 ) + S~(p k+1/2 ,q k+1 ) (56) 

h 

= Pk+l/2 Qk+i ~ Pk+l/2 Qk~ 7; { H (q k ,p k+1 / 2 ) + H(q k+ i,p k+1/2 )) , 
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and a extremal condition in the intermediate variable Pk+i/2- That is, 

OS , OS , 

d«S = -7T- dq k + - dg fc+ i 

dq k oq k+ i 

\ <9g fc dp k+1/2 J dq k J fc 
+ I — ^ — I — — — f d Qk+i, 



dq k+ i \dp k+ i/ 2 dp k+1/2 J dq k+1 
which leads to 

iw^ 11 (57) 

(put in another way, S(q k , q k +i) is not a function of p k +\/2 and consequently 
its partial derivative with respect to this variable should vanish). In other 
words, we obtain the first equation in (|5ip by —p k = ^is^dh+ll ; the third 

one by Pk+l = 85( ^ fc + l) , and the second one by a5 ^ gfc + l) = 0. 

As shown in Theorem 15.41 S(q k ,q k+ i) is projectable onto L d , while con- 
dition d£7J) provides C d C M x M. 



6.3.1. Regular systems. Consider the usual mechanical Hamiltonian 

H(q,p) = ±pM- 1 p T + V(q), 



where M is a symmetric regular n x n matrix. From the second equation 

-1/2 



in (f5"Tj) is easy to check that P^ +l / 2 = M( Qk+1 h Qk ). Hence, projecting ([56 



onto M x M we arrive to 

In the expression just above, we clearly recognize the discretization 

T , \ h , q k+ \-q k h Qk+i~qk^ 
L d (q k ,q k+1 ) = -L(q k , ) + -L(q k+1 , 



2 w " h ' 2 h 
i 

2 



for the usual mechanical Lagrangian L(g, g) = \ q T M q — V(q). 



6.3.2. Martinet structure. Consider again the sub-Riemannian Martinet struc- 
ture in (|4"2"|h Recall that q = (x, y, z) T G M 3 and p = (p x ,Py,Pz) € M 3 . From 
([57|) we obtain 

a; fc+ i = a; fe + ^ ( ('(p fc+1 / 2 )a + (Pk+i/2)z^f) + f(pfc+i/ 2 )» + (Pfc+i/2) z ^f i 



/ (Pfc+l/2)y (Pfc+i/a)* 
2/fc+i = y*: + - < 77—7, — vi + 



2 {(l + f3x k ) 2 (l + /3x k+1 )\ 

Zk+1 = 2 fc + ^ j^(p fe+ i/2)* + (pfc + i/ 2 ) z y^ y + + (Pfe+i^)^^-^ ^y 1 } 

As above, we consider a too steps (of /i/2 size) interpretation of this sub- 
Riemannian system through equations (|52p and (|53p in the following manner 
(we denote p k +i/2 by p for sake of simplicity): 
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(qk,Pk) -> (<?fc+l/2,Pfc+l/2 



Xk+1/2 = Xk + | (p X + Pz\ 



Vk+i/2 = Vk + h { {1+ % k) 2 } 



Zk+l/2 = Zk + ^ (px+Pz^-) \ 



I z k+\/2~ z k \ _ Vk ( x k+\/2~ x k 
• \ h/2 ) ~ 2 \ h/2 



Ldilk, Qk+l/2) 



(Qk+l/2>Pk+l/2) (Qk+l,Pk+l) 



Xk+l = Xk+1/2 + I (px +Pz^Y 1 ^j 



2 + h { {l+P X +l) 2 } 



Vk+i = Vk+i/ 



4 



„ , z k+1 -z k+1/2 \ _ y k+1 f x k+1 -x k+1/2 

■ [ h / 2 ) — 2 \ h/2 

-^d (ftc+1/2, Qk+l) 



As shown in the table above, both discrete submanifolds are respectively 
fined by the constraints (*±jgp) = # (^^) and (^2^) = 

. Taking into account equations flU]) and ([55]) in § 16.31 



2 ^ h/2 

define the two discrete Lagrangian functions (implicitly expressed) as 

L d (?fe) 9fc+l/2) = Pfc+1/2 <37c+i/2 ~Pk+i/2Qk - -^H(q k ,p k+1 / 2 ), 

with qk+1/2 = Qk + ^H p (q k ,p k+1 / 2 ) (recall that this expression corresponds 
to the generating function after projecting onto M x M). On the other hand 

L d 0?fc+i/2><?fc+i) = - Pfc+l/2 9jfc+l/2 + Pfc+1/2 <7fc+i - -^H{qk+i,Pk+i/2)i 

with g fc +i = qk+1/2 + ^H p (q k+ i,p k+l / 2 ). 

From the previous expressions and following the discussion in [30] § 2.5.1, 
we divide each step (q k ,q k+ \) into 2 substeps (qk = q k , q\ = qk+1/2) an d 
(% = Qk+i/2 j <7fc+i = Qfc+i)- Take the discrete action sum 

JV-l 

©.({^ir 1 ) = E( L M>^ + ^(«+i)) 

A:=0 
JV 

= ^2 ( L d (Qk, qk+1/2) + L d (q k+1/2 ,qk+i)) ■ 

k=0 

The corresponding Euler-Lagrange equations, resulting from requiring this 
action to be stationary, pair both neighbouring discrete Lagrangians to- 
gether to give 

D 2 L+(q° k .ql) + D 1 L^,^) = 0, 
D 2 L-(qlq 2 k ) + D 1 L+(q° k+1 , q\ +1 ) = 0. 
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The equations just above completely determine the discrete dynamics in 
equations (|52j) and (|53p for the Martinet sub-Riemannian system (|42p . In 
other words, the discrete scheme in ([52]) and ([53]) could be understood as 
two discrete Hamiltonian flows F L + and F L -, each one of half-step h/2, 

d d 

respectively generated by the generating functions and L~^. The map 
over the entire time-step h is thus the composition of the maps F L - o F L + . 



7. Conclusions and future work 

In this paper, we have carefully studied the relationship between Hamil- 
tonian dynamics and constrained variational calculus, showing that under 
natural regularity conditions both are equivalent. The analysis is easily 
extended to the discrete case, allowing us to explore the applications to 
the construction of symplectic integrators. Moreover, we have analyzed in 
parallel the case of classical nonholonomic mechanics in the discrete and con- 
tinuous cases. Our technique allows us to extend the comparison algorithms 
between nonholonomic and constrained variational cases to the discrete pic- 
ture, finding additionally a new and simpler point of view in the continuous 
case. 

Some interesting lines that we want to explore in future research are the 
following ones. First, many higher-order symplectic methods are obtained 
using composition of methods of lower-order. In our approach, this notion 
seems to be related with the notion of composability of canonical relations 
(which may fail even to be a manifold without appropriate transversality 
conditions |37|). This study is a promising line for numerical simulation of 
Hamiltonian dynamics and also constrained systems. 

One of the main problems of discrete nonholonomic mechanics is the lack 
of preservation of geometric structures (non-preservation of the nonholo- 
nomic bracket, non-preservation of the nonholonomic momentum in gen- 
eral) mimicking the non-preservation of the continuous nonholonomic sys- 
tem. Then, it is difficult to compare from this geometric perspective the 
nonholonomic integrators obtained from discrete nonholonomic mechanics 
with standard methods. One possibility is to use our comparison algorithms 
detecting if one particular nonholonomic integrator is preserving the com- 
mon solutions for the continuous nonholonomic problem and its associated 
constrained variational problem. 

Other interesting case that we want to explore is the extension of our 
theory to reduced systems using the geometric framework given by the Lie 
algebroid and Lie groupoid formalisms [28j . For instance, in the discrete 
case, our intention is to derive the dynamics using also Lagrangian sub- 
manifolds of the so-called tangent groupoid [11]. In this context we will 
start with a submanifold N of a Lie groupoid G M, and a discrete La- 
grangian : N — y R. Following Theorem 12. 11 we introduce the Lagrangian 
submanifold of the cotangent groupoid T*G. The Hamiltonian side 
is determined by a Poisson flow defined on A*G, the dual bundle of the 
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associated Lie algebroid AG to G. Therefore, we have the following scheme: 




which, in principle, will allow us to extend the theory presented in this 
paper. 
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